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Abstract. We study how a chiral tricritical point (TCP) on QCD phase diagram is affected by the 
imbalance of up and down quark densities (isospin density), using the generalized Ginzburg-Landau 
Q-i (GL) approach. The resulting phase diagram near TCP shows a rich fine structure which includes 

inhomogeneities of both the chiral and the charged pion condensations. It turns out that the TCP 
splits into multicritical points. 
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INTRODUCTION 

43 ■ 

Study of the phase structures of QCD at finite temperature and quark density is intriguing 
subject. In particular, it is suggested there exists a critical end point (CEP) in the 2- 
flavor QCD phase diagram [1], which leads to a tricritical point (TCP) in the chiral limit, 
although the location of the CEP is not well determined. According to recent work based 
on the Ginzburg-Landau (GL) approach [2, 3] and effective models [4] for the isospin 
symmetric 2-flavor quark matter, the TCP is shown to be replaced with the Lifshitz point, 
at which three different phases meet; the Wigner phase, the chiral symmetry broken 
phase with the spatially homogeneous chiral condensate (HCC), and the broken phase 
with the inhomogeneous chiral condensate (ICC). 

How the phase structures near the TCP will change under realistic situations for 
the hadronic matter? We focus on the effect of the isospin asymmetry, which may be 
realized in the quark matter under the charge neutrality condition or in the realistic 
heavy-ion collisions. There are only a few model-based works on the TCP with the 
isospin asymmetry [5]. On the other hand, the charged pion condensed phase is known 
to show up as the true ground state of QCD at large isospin density [6], We therefore 
try to present a systematic model-independent approach based on the generalized GL 
method to study effects of the finite isospin density on the TCP and its neighborhood. 

In this report we focus on the phase structure in the vicinity of TCP for the massless u 
and d quark matter based on our recent work [7]. Taking into account isospin dependent 
contributions to the GL functional, we obtain new fine structures near TCP. As a result, 
the charged pion condensates extend over sizable domains in the GL parameter space. 
The finite isospin density also brings about a splitting of TCP into four independent 
critical points, which will be shown below. 



SYMMETRY BREAKING DUE TO THE ISOSPIN DENSITY 



In the vicinity of the chiral TCP, the local free energy for the chiral phase transition can 
be expressed in terms of the chiral 4-vector (j) = (<J,7t l ) as order parameters, where o 
is proportional to the chiral condensate (qq) and K l the pion condensates (qi^T'q) (i = 
1, 2, 3) with doublet quark fields q = (w, d). In order to provide a minimal description of 
TCP, we expand the GL potential up to the sixth order of the order parameters, and their 
spatial derivatives as follows, 

Q GL [a(x)y(x)] = f 2 + f (0 2 ) 2 + ^(V0) 2 +|(0 2 )3 + ^(0,V0) 2 



[0 2 (V0) 2 - (0, V0) 2 ] + -P(A0) 2 + 8w c + 8(0! , 
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where 0^/3, are the GL parameters, 8co c stands for the contribution from current quark 
mass m c , and 8(Oj expresses the contribution from the finite isospin density. Eq. (1) 
except 8(O c and 8(Oj is invariant under the chiral SU(2) L x SU(2) W transformation. The 
explicit chiral symmetry breaking term, 8(O c , may be given by —ho with a small constant 
h. We set h = throughout this report. 

For the isospin symmetry breaking term, 8(D/, we consider the GL potential written 
with the charged pion condensates, = (ft 1 , 7T 2 ), which is given by 

8co, = |(0 2 + 1(0 4 + % if ~ (O 2 ] W) 2 + %(V<) 2 • (2) 

Here, we retain the 7t£ contributions in addition to the leading 7T 2 terms, because we 
apply our GL analysis to the region where the order parameters become comparable with 
isospin density. Without loss of generality, we set 7t l = 7t, n 2 = n 3 = from symmetry. 

The GL potential Eq. (1), (2) contains parameters a ; ,/3,, which are, in principle, the 
function of temperature and quark/isospin chemical potentials. In order to reduce the 
number of the independent parameters in Eq. (1), (2), we make use of the feedback of the 
quark loops to the potential energy, assuming that they are the dominant contributions. 
We systematically relate the GL parameters with the n-th quark loop integrals, which 
certainly depend on the u and d quark chemical potential, ji u and ji^. Expanding the 
resulting integrals up to the 2nd order of the perturbative isospin chemical potential 
111 = {p, u — Hd)/2, we can show the GL parameters, a(/i/), /3(/l/) are proportional to 
a(f±i = 0) or ju 2 a(jli = 0) [7]. Hence, we finally obtain the potential: 

Ogl(jU/) = «f (J 2 /2 + (af> - ^af ] /2)K 2 /2 

+(4 0) +ti4 0) ) [tf 4 + aV+(a') 2 ]/4 + af [^ 4 + (7 2 ^ 2 + (^) 2 ] /4 (3) 
+a^ 0) [(a 2 + k 2 ) 3 + 5(<jo' + kk') 2 + 3(<jk' - kg') 2 + [(a") 2 + (k") 2 }/2] /6 , 

where = = 0) . We restrict the analysis to one-dimensional structures and 
accordingly a primed quantity in the above formula means its derivative with respect 

to z-direction, for instance, a' = d z a. We simply choose = 1 to set the energy 
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FIGURE 1. Phase structure near the TCP for /X/ ^ 0, h = 0. The solid and dashed curves indicate 
boundaries of the first and second order phase transitions, respectively, (a) Phase diagram with the HCC 
and HPC: Asymptotic behaviors of the first order critical line separating the HCC phase from the HPC 
phase is analytically derived: Cfe -» 3jlt/a4/4 + ju//24 + ff{jxf) as U\ -> while a,2 -> nfo^/l + a^/S 
as 0:4 — > 0. (b) Phase diagram with the ICC and IPC in addition to homogeneous condensates: The second 
order critical line connecting TCP' to Q is characterized by (X2 — 3 ( OC4 + (if ) 2 /8 . The second order critical 
line linking TCP' to R is given by a-i = 5(0:4 + M/ 2 ) 2 /36. The first order critical line between the ICC and 
IPC phases (a curve QR) is obtained by a numerical computation. 



TABLE 1. The location of chiral multicritical points; 
TCP' and P in Fig. 1(a), Q, R in Fig. 1(b). 
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scale. By virtue of the appropriate scaling of the GL parameters [3], the potential 
Eq. (3) is characterized by two independent parameters, {0C2/CC4, /x/ / 1 CC4 1 } . Hence, we 

shall describe the phase diagram in the (a^\a^) parameter space with the given jij. 
Hereafter, we suppress the subscript (0) for as to avoid notational confusion. 



POSSIBLE PHASE STRUCTURE AT FINITE ISOSPIN DENSITY 

Let us first assume the spatially constant distribution for the chiral and charged pion 
condensates. We show the phase diagram near the TCP in Fig. 1 (a), and locations of 
multicritical points on the GL parameter space in TABLE 1 . In this case, the ground 
state favors either HCC phase or homogeneous pion condensed (HPC) phase. There is 
no room for the coexistence phase of % and o. Compared with the result for the /i/ = 
matter, large domains of the phase diagram is occupied by the HPC phase. This drastic 
change is due to the presence of (0C2 — jU 7 2 0(4/2) 7T 2 /2 term in the GL potential Eq. (3), 



which favors the HPC phase for arbitrary small /!/ with CC4 > 0. In addition, while the 
TCP is located at the origin of the phase diagram for the /i/ = matter, the finite /!/ 
brings about a shift of the location of TCP to (0, -/i/ 2 ), newly labeled by TCP'. Origin 
of this shift is easily understood by looking at the coefficient of the o 4 term in Eq.(3), 

a 4°" > + M/ 2 &6> (ot6 = 1). Our result is consistent with the model-based calculations [5]. 

Next, we consider possible inhomogeneous structures for the chiral and charged pion 
condensates with the finite isospin chemical potential. Following the previous work, 
we consider the one-dimensional solitonic structure expressed by the Jacobi's elliptic 
function "sn", for the inhomogeneous phases of o and n condensates [3]. In fact, the 
Jacobi's elliptic function is shown to be a solution of the Euler-Lagrange equation for 
the GL potential in the chiral limit. We assume c(x) = ^Jvkvc\(kz\ v), K = for the ICC 
phase, and o = 0, n(x) = v /vfcsn(fe;v) for inhomogeneous charged pion condensed 
(IPC) phase, where the positive parameters V and k are determined by the variational 
principle in each phase. The resulting phase diagram is depicted in Fig. 1 (b). 

The TCP is now replaced with the Lifshitz tricritical point. Compared with isospin 
symmetric matter, a large part of the ICC phase is taken over by the IPC. This is 
because the spatial derivative term of (n') 2 gets lower energy compared with the (a') 2 
contribution as clearly seen in Eq. (3). The boundary between ICC and IPC also provides 
two new critical points, Q and R, listed in Table 1. 

For completeness, we also consider possible o-K coexistence phases. We take the 
ansatz of chiral kink spiral: c(x) = mcn(fe; v), n(x) = msn(fe;v) with m a new vari- 
ational parameter. This includes the standard chiral spiral c(x) = mcos (fe), 7r(x) = 
msin(fe) [8] in a particular limit v — > 0. We always find either IPC or ICC being fa- 
vored, indicating that the a-lt coexistence state may be ruled out near TCP. 

To summarize, we studied the phase structures in the proximity of TCP at finite /!/, 
using the model-independent GL approach. The sizable domain of the phase diagram 
was found to be occupied by the charged pion condensation. We found not only the 
shift of the chiral TCP, but also the appearance of the new critical points. The new fine 
structure obtained in the GL parameter space may be mapped onto the low temperature 
and intermediate density of the (T,/i) QCD phase diagram, as suggested in the previ- 
ous work [4]. The inhomogeneous charged pion condensate of the quark matter may 
smoothly continue to the classical pion condensate discussed in the nuclear matter [9]. 
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